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Abstract 

By making use of a non-perturbative definition of a momentum diffusion coefficient as well 
as Heavy Meson Chiral Perturbation Theory, we investigate the Brownian motion and ki- 
netic equilibration of heavy quark flavours deep in the confined phase. It appears that the 
momentum diffusion coefficient can be expressed in terms of known low-energy constants; it 
increases rapidly at temperatures above 50 MeV, behaving as ~ T 1 /F% for ^ « T « F w , 
where m n and F n are the pion mass and decay constant, respectively. The early increase may 
suggest a broad peak in n/T 3 around the QCD crossover. For a more detailed understanding 
the computation could be generalized in a number of ways. 
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1. Introduction 



An ideal scenario for what could happen in a heavy ion collision is that light quarks and 
gluons form a rapidly expanding thermalized medium, and that there are some "probes" 
available, whose properties are affected by the medium in a significant yet tractable way. 
Among the most attractive probe candidates are heavy quarks (charm and bottom quarks 
as well as their antiparticles) , which can be copiously produced in an initial hard process. 
After a while the heavy quarks decay, but modifications on their behaviour caused by the 
thermal medium could conceivably be deduced from the experimentally observed transverse 
momentum distributions and azimuthal anisotropics of the leptonic decay products [U [2] . 

With such motivations in mind, a significant body of work has been carried out during the 
last 20 years or so, concerning the effects that a thermal medium can have on the propagation 
of heavy quarks J3[ d] . Roughly, the initial stages are characterized by radiative energy loss 
(bremsstrahlung), which may slow down the heavy quarks; the final stages are characterized 
by elastic scatterings, which cause collisional energy loss but also produce random kicks 
corresponding to Brownian motion. Assuming, idealistically, that the thermal system is 
spatially large enough for all of these processes to take place, they have been described by a 
number of related physical concepts and observables, such as energy loss (dE/dx), stopping 
distance, jet quenching, (momentum) diffusion, drag, or kinetic equilibration. 

In the present paper, we focus on a late stage of the above scenario, in which the heavy 
quarks are practically at rest with respect to the thermal medium but also undergo Brown- 
ian motion. (Possible extensions to more general kinematics are outlined in the conclusions.) 
More specifically, we are interested in an observable called the momentum diffusion coeffi- 
cient, k, which characterizes the random force acting on the heavy quarks and, through a 
fluctuation-dissipation relation, also determines their kinetic equilibration rate. 

Recently, a number of theoretical works have addressed the same problem. For instance, 
k has been computed to leading [5] and next-to-leading [6] order in the weak-coupling ex- 
pansion. It has also been given a non-perturbative definition [7] in the framework of Heavy 
Quark Effective Theory [8]-[12]; this was partly inspired by computations through AdS/CFT 
techniques in the large- N c limit of strongly coupled M = 4 Super- Yang-Mills theory, which 
suggest a larger k than in leading-order QCD [131 [HI [15] . (Similar conclusions have also been 
reached for AdS models resembling QCD; see e.g. ref. [16] and references therein.) Numerical 
simulations have been carried out within so-called classical lattice gauge theory, suggesting 
again a value larger than indicated by the weak-coupling expansion [T7] . All of this makes a 
strong case for attacking the problem with lattice simulations, a challenge that may be less 
daunting that the determination of many other "transport coefficients", such as viscosities, 
because of the simpler structure of the pertinent spectral function [18]; indeed the very first 
numerical attempts look rather promising [19J . The inherent uncertainties related to analytic 
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continuation from Euclidean numerical data have also been looked into, and it appears that 
it might be feasible to carry through the program at least on the qualitative level [20]. Once 
a value of k is available, it can be incorporated in hydro dynamical simulations to yield results 
relevant for experiment (cf. e.g. ref. [5]); there is a growing body of such works under way. 

The purpose of the present paper is to make use of the non-perturbative definition of k 
introduced in ref. [7], and to evaluate it deep in the confined phase. This is possible in QCD 
because, due to chiral symmetry breaking^ the infrared dynamics of the confined phase can 
be parametrized with a small number of "low-energy constants". In fact, in a certain limit, 
we find that k is fully determined in terms of the pion decay constant and the pion mass. 
There is, of course, a long history to applying Chiral Perturbation Theory (not to mention 
hadronic models) for the computation of various thermodynamic properties of the pion gas, 
see e.g. ref. |21j ; past developments and the non-trivial theoretical challenges that are related 
particularly to transport coefficients have recently been discussed in the context of the bulk 
viscosity of strongly interacting matter in ref. |22j . 

The paper is organized as follows. The non-perturbative definition of the momentum 
diffusion coefficient is reviewed in sec. [2j the chiral effective theory relevant for handling 
heavy-light mesons in the confined phase is described in sec. [3j and the computation of the 
momentum diffusion coefficient is presented in sec. [U Some conclusions comprise sec. 

2. Basic physics of momentum diffusion 

Considering time scales short compared with the life-time of the heavy quarks (or heavy-light 
mesons) but long compared with those of microscopic processes involving gluons and light 
quarks (or light mesons), which have a momentum p ~ T ~ 100 MeV, where T denotes the 
temperature, the dynamics of the heavy degrees of freedom, with a mass M » T, presumably 
resembles Brownian motion. If so, it can be described by the Langevin equation, with the 
role of the stochastic noise being played by a force induced by QCD-mediated collisions with 
gluons and light quarks (or light mesons), which are present with an abundant number density 
n ~ T 3 . The corresponding equations of motion have the form 

Pfc(t) = -VoPkit) +&(t) , (2.1) 
«&(t)&(f)» = K6 kl 6(t-t'), «£ fe (i)» = 0, (2.2) 

where pt is the momentum of the heavy objects (k = 1,2,3); ^ is a Gaussian stochastic 
noise; and ((•••)) denotes an average over the noise. According to eq. (|2.2p the momentum 

x The number of light flavours is assumed non-zero, N{ > 0; it is unclear whether k can be given a sensible 
meaning in the confined phase for Nf = [7] even if the related Euclidean correlator appears to exist [18] . 
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diffusion coefficient, k, characterizes the auto-correlator of the force, 

1 f 00 

« = o dt j]«e*(*)efc(o)», (2.3) 

a J-co fc 

whereas the coefficient rj D appearing in eq. (|2.ip is referred to as the "kinetic equilibration 
rate" or the "drag coefficient". As is well-known, in classical statistical physics the two 
can be fluctuation-dissipation related to each other: rj D ~ n/(2TM), where we interpret M 
as a specific mass definition (more precisely the heavy quark "kinetic" mass; we assume a 
regularization scheme respecting Lorentz symmetry, and then M is equal to the "rest" or 
"pole" mass; it should be fixed non-perturbatively). 

Now, it was argued in ref. |7] that the classical definition of k as an autocorrelator of forces, 
eq. (|2.3p . can "naturally" be extended to QCD. Suppose that we know the heavy quark Hamil- 
tonian, H, as well as the Noether current associated with the U(l) flavour symmetry, . 
More specifically, H and J® are needed up to 0(M ) in an expansion in a large M, whereas 
the J are needed up to 0(l/M). With these operators we can define a "susceptibility" 
related to the conserved charge, 



J (j%x)J%0)) T , = 1, J=Jd 3 x, (2.4) 



and the "acceleration" associated with the spatial components, 



[H,J k ]+^~, (2.5) 



where the partial derivative acts on possible background fields. Consequently. 

'1 (dj h (t,x) dj k (t',0) 



f V limf lim ^ fdte^-^f (- 
3 o[m^x 00 J-oo Jx\2 



K 

■A 

fe=l 



dt ' At' 



(2.6) 



Since ML dS k /dt represents, according to Newton's law, a force, this definition is indeed 
a generalization of eq. (j2 . 3|) . The ordering of the various limits requires, however, a careful 
analysis [7]; the situation simplifies only if a dependence dj k /dt ~ 1/M can be factored out 
and cancelled. We also note that, as dictated by standard relations between various time 
orderings at finite temperatures [23J, 

k = lim , (2 . 7) 

where the spectral function p E is defined as in eq. (|2.6p but with a commutator replacing the 
anticommutator. (The notation (. . .) T refers to the usual thermal average.) 
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3. Heavy Meson Chiral Perturbation Theory 

The goal now would be to evaluate eq. (|2.6[) at low temperatures. This can be achieved by 
making use of an effective field theory that is valid in the regime considered and allows us 
to define the operators entering the definition. Given that the framework may be unfamiliar 
within the finite-temperature community, we give a few ingredients in the following, although 
no attempt is made at a comprehensive review (see refs. \24\ 1251 [26] for introductions). 

The "usual" chiral Lagrangian describing the pseudo Nambu-Goldstone bosons of chiral 
symmetry breaking has the form [27] 

F 2 Y 
£ xPT = _Tr (dHJdpUl) + — Tr (M^U + U^M) + ... , (3.1) 

where F is the pion decay constant in the chiral limit (naively F ~ 93 MeV), E is the 
chiral condensate, Ai is the quark mass matrix, and U e SU(iVf) is the Goldstone field. For 
simplicity we take the mass matrix to be of the form M. = m 1 in the following, with m s R; 
then the pion mass squared is m 2 = 2mS /F 2 + 0{l/F A ). Parametrizing U = exp(^), the 
tree-level propagator is 

(UW tcM > (0) = \ {KAc ~ r^5 ab 5 cd )A(x - y; m 2 n ) , (3.2) 

where A is a scalar propagator. It is important to keep in mind that for the A4 chosen, 
eq. (|3.ip contains no terms cubic in £, i.e. no three-pion interactions. (Interactions among 
odd numbers of pions are contained in the Wess-Zumino-Witten term, but it is suppressed 
by a large number of derivatives.) 

The next step is to supplement the chiral Lagrangian with a piece describing heavy-light 
mesons. The resulting effective description is called Heavy Meson Chiral Perturbation Theory 
(HMxPT) [281 ESI EQ] . 

To start with we note that, unlike Heavy Quark Effective Theory (HQET) [5]— [12j. which 
describes all states containing a single heavy quark or antiquark (provided that the associated 
gluons and light quarks are soft), the mesonic HM%PT only contains a specific subset of such 
states, namely parity-odd pseudoscalar and vector mesons (D and D* for charm, B and B* 
for bottom), with total spin or 1. The standard convention is to denote the scalar by P a 
and the vector by Pa^, where o is a light flavour index (a = 1, . . . , JVf); to streamline the 
notation we define Qa = Pa 11 in the following. Although Qa is written as a four-vector, it 
only contains three independent components, cf. eq. (|3.10l) below. 

The form of HM^PT is dictated by symmetries, of which there are many. First of all there 
is a U(l) symmetry related to the conserved heavy quark number; this yields the Noether 
current alluded to above. Second, HQET at O(M ) displays a heavy quark spin symmetry, 
since Pauli matrices first appear at 0(1/M); this implies a relation between the fields P a and 
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Qa (cf. eq. (|4.2p below). Third, the light-quark index a enjoys a specific transformation 
property under the SU(iVf)£xSU(./Vf)^ chiral symmetry. Fourth, even though the effective 
Lagrangian is non-relativistic, its origin in relativistic QCD implies that it also remembers 
something about proper Lorentz symmetry, provided that this has not been broken through 
regularization. This is usually implemented by introducing an arbitrary parameter, a four- 
velocity v^, and an associated "velocity reparametrization invariance" |31j . Finally, there are 
the usual discrete C, P, and T symmetries of QCD. 

It turns out to be non-trivial to implement all the symmetries in a convenient way. In 
fact, in order to define a simple parity transformation it is not practical to use the Goldstone 
field U in the part of the Lagrangian involving the heavy mesons, but rather a "coset field" 
(General formal considerations can be found in refs. |33|. I34j.) If U transforms in 
SU(iV f ) L xSU(iV f ) fl as 

U->LUltf, (3.3) 
then the field \/U can be assigned the transformation 

VU -> LVUW^ and VU -» WVUrf , (3.4) 

where the complicated (space-time dependent) field W is defined through eqs. (|3.4p . We can 
subsequently introduce the traceless and Hermitean fields 




which transform as 

-> WV„W^ + i WdpWt , A,, -» WA„W^ . (3.7) 
Note that in the chiral expansion, inserting \f\J = exp(^), we get 




i.e. V M couples to an even number of pions and A^ to an odd number. 

The heavy meson fields P a and Qa are normally assembled into a 4 x 4 matrix H a and its 
conjugate H a = j°hIj°: 

H a = l±t (q a + i 75 p a ) , s a = (gl + i7sP t) 1±1 , (3.9) 

where, following a frequent convention, (...y denotes complex conjugation acting on -fa; Qa- 
Here is a four- velocity with v ■ v = 1. The vector meson field is constrained by 

v ■ Q a = , (3.10) 

2 We follow the notation of e.g. ref. [32], where HMxPT was used in connection with a lattice investigation. 
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which implies for instance that tfQ a = —Qai- The chiral transformation can be assigned as 

H->HW*, H^WH. (3.11) 

Defining 

the effective Lagrangian at 0(M ) has the form 

4m X PT = -» Tr (#« " • #&) + M p TV (# a # a ) + ^ Tr {H a H b 4 ba ^) 
+ ai Tr (H a H b )(VUM^VU + VU^MVU^) ba 

+ a[Tv(H a H a )(VUM^VU + Vu t MVU t ) bb + ... , (3.13) 

where the trace is over Dirac matrices, and ". . ." denotes terms of higher order in the chiral 
expansion. The coefficient is a new dimensionless low-energy constant, estimated from 
e.g. D* — > Dtt decays [26J or lattice determinations (e.g. refs. [351 [36]) to be around g n ~ 0.5. 

The mass term on the first row of eq. (|3.13p . proportional to M p , is often not shown, 
because it can be shifted away by a time-dependent phase transformation. However, for the 
finite-temperature application that will be discussed in the next section, it is useful to keep 
M p explicit. The terms proportional to the low-energy constants cri,a[ play a role, at leading 
order, for the mass spectrum: a[ shifts M p by a flavour-independent amount, whereas a% 
would break the flavour symmetry if A4 were non-degenerate. In our case both of these effects 
can be accounted for through M p , so we omit <J\,o\ in the following. The spin symmetry 
between P a and Qa is only broken at 0(1/M); although we presently turn to such effects, 
the "trivial" mass splitting that they induce will be considered to already be accounted for 
by the present remarks, and we denote the vector meson mass by Mq. 

Now, we turn to terms of 0(l/M). Some of these can be deduced from eq. f)3. 13j) through 
reparametrization invariance [31] . This is a transformation which leaves invariant the "origi- 
nal" momentum operator = Mv^ + W^ as well as scalar products such as v-v = 1. Writing 
Vfi — > Vfj, + e M /M with v • e = 0, which implies iD^ — > iD^ — e^, we observe in particular that 
the combination 

O ba = ivV ba --^(V-V) ba (3.14) 

is invariant only for a± = 1. This is clearly a reflection of Lorentz symmetry, and means that 
the coefficient of a particular 0(l/M)-operator is fixed in terms of that of iv ■ T> ba at O(M ). 
(Whether it is M or M p that appears in the denominator of eq. (|3. 14[) plays no role, because 
their difference amounts to an effect of 0(1/M 2 ).) So, 

4m x pt = ^ Tr & ■ V ^a H b ] + ... . (3.15) 

The term proportional to g n in eq. (|3.13p also implies the existence of terms of 0(g n /M) with 
specific coefficients, but the argument is not quite as transparent as that with eq. (|3.14j) : we 
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postpone a discussion to after eq. (|4.2p . Analyses of all terms of 0(1/M), also those not fixed 
by reparametrization invariance, have been carried out in refs. |3T|, I38j. 



4. Main computation 



We now apply the theory defined by eqs. (|3.13p . (|3.15p to estimate the momentum diffusion 
coefficient as defined by eq. (12.6p . The computation carried out is rather straightforward and 
simple-minded, but the result is interesting so we present the steps in some detail. 



4.1. Force operator 



We start by rewriting the sum of eqs. f|3. 13|) . f|3. 15|) after carrying out the Dirac traces, setting 
v — > (1, 0), and adding terms of 0(g n /M) through an argument to be discussed presently (the 

indices k,l,m are spatial; indices are raised and lowered with the metric (h ); repeated 

indices are summed over; and ". . ." indicates terms of higher order in the chiral expansion): 



-HMxPT 



M P PjlP„ 



2[Ql k iV° ba Q bk -M Q Ql k Q ak 



^[p!(V l V l ) ba P b + Ql k (V l V l ) ba Q bk ] 

j [2PlV k cb Q ck + 2Ql i D k ac P b + e klm (Ql k V l cb Q cm - Q^LcQhn 



M 
9_ 

M 



A? 



ba 



(4.1) 



In the part of 0(1/M) only terms with spatial derivatives acting on P and Q have been 
shown, because these are the only ones contributing to J k . Also, QVd ac = (V^Q^. 
The heavy quark spin symmetry amounts to the transformations 



SP„ 



-a t Q 



al i 



SQ, 



ak 



a k Pa + e klm a l Qa 



a k eR, k,l,m e {1,2,3} . (4.2) 



Terms of O(M ) in eq. (|4.ip are easily seen to be invariant under this transformation; the 
mass terms are invariant only if Mq = M p , i.e. the difference Mq — M p must be of 0(1/M), 
as mentioned above. Even though this is not the case in general, the terms of 0(1/ M) shown 
explicitly in eq. ()4.ip are also invariant under eq. f)4.2|) (the term on the fourth row is invariant 
only up to total derivatives and terms containing covariant derivatives acting on A ba , which 
have been omitted because they do not contribute to J k \ The invariance exists because 
these terms are related to invariant terms of O(M ). 

The fourth row of eq. (14. ip . with a coefficient g n /M, is new with respect to eq. (|3.15p . 
It is connected to the second row of eq. (14. 1|) through reparametrization invariance; to be 
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Figure 1: Some of the scatterings experienced by P- mesons almost at rest with respect to a pion gas 
(P and P* could represent B and B*-mesons, respectively). The process (a) is kinematically allowed 
for D-mesons, but not for £?-mesons. There is no i-channel scattering on thermal pions, because the 
Lagrangian in cq. (|3.ip contains no three-pion vertex and that in eq. (|4.1[) contains no PP7r-vertex. 
Electromagnetic interactions, such as scatterings on blackbody photons, have been omitted. 



concrete, one could imagine that the vector field of an original relativistic theory, let us 
denote it with Q, satisfy a transversality constraint of the type P^Q^ = rather than 
eq. (|3.10p . So, recalling that P^ = Mv^ + iD^, every appearance of can be replaced with 
Qt 1 = QA* _ ■ Q/M, where Q» does satisfy eq. (I3T01) [31J. This yields terms of 0(g- K /M) 
with fixed coefficients. The 0(1/ M) structures shown in eq. (|4. 1[) that have been determined 
through reparametrization invariance are the only ones from the full list [38] that contain 
spatial derivatives acting on P a , Q a u an d therefore contribute to J k . (An explicit crosscheck 
on the logic presented could be obtained by starting from a genuinely relativistic formulation, 
such as the ones in refs. [37l [39] , and taking the non-relativistic limit only in the end, because 
then all consequences of Lorentz symmetry are automatically respected. In the following the 
terms of 0(g n /M) are included for illustration but are omitted from the final result.) 

Now, all terms proportional to g w couple to A^ a and thus to an odd number of pions, 
whereas terms without g w couple to an even number of pions (cf. eq. (|3.8|) ). Also, g n nec- 
essarily mixes the fields P and Q\ or Q and Q^. This makes the terms proportional to g n 
in general the most important ones from the point of view of zero-temperature phenomenol- 
ogy, such as the decays of the D* and B* mesons [26]; a relevant amplitude is illustrated 
in fig. [D^a). However, if all particles are on-shell, this process is kinematically forbidden for 
S-mesons, because the mass difference between B and B*, scaling as 1/M, is below the pion 
mass. Therefore the dominant processes for U-mesons are of the types shown in figs. QJb) 
andCDlc), and the scattering amplitude is 0(1/F 2 ). 

Next, we obtain from eq. (|4. 11) the U(l) Noether current: 

J° = 2PlP a + (P a ^Q ak )+o(^) , (4.3) 

^ = ^ P a dkP a-S k PlPa)-^PlV k ba P b ]+(P a ^Q al ) 

+ 2 -Jf[ P lQbk - Ql k P b + £lkmQilQ b m]A° b a + O Q^j . (4.4) 
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In order to determine the acceleration from eq. (|2.5p . we would need the Hamiltonian; how- 
ever, equivalent information should be contained in classical equations of motion, which read 

id P a = M P P a + VlP b -ig,A l ba Q bl + o(^) , (4.5) 
id Q ak = M Q Q ak + Vl Q bk + ig n (A k ba P b - e klm A l ba Q bm ) + O (— ) . (4.6) 

Note that when acting on J~°, J k , the terms of 0(M P , Mq) from here cancel out (apart 
from a possible remainder proportional to Mq — M p ~ 1/M), so that effectively the time 
derivatives count as terms of O(M ). 

Combining eqs. (|4.4p - (|4,6p . and noting that in eq. (|2,6p a spatial average can be taken over 
the currents and that therefore partial integrations are allowed with respect to d k , we finally 
obtain <9° $ x i7 fc . Terms do proliferate quite a bit: defining 

iF% = [^ T ,V^ b = i(d^ a -d^ a + iV^ c -iV^), (4.7) 

= [V» T ,iA» T ] ab = l (d»Al + l V? a Al-iA»X c ) , (4-8) 
%H% = [iA» T ,iA» T ] ab = i(iAZ a A» bc -iA» a A£ c ) , (4.9) 
iH% = {iA^ T M vT U=i(iA^Al c + iA v ca A^, (4.10) 

where {..) T denotes a transpose with respect to flavour indices, we obtain the structure 

d°[j k = -| j\pHF k0 - £H M )P + Q}(F k0 - glH^Q, 
7tt [p^ G lk Qi — Q\ G lk P + e mln Q^ m G lk Q n 

ln [ P t (G oo _ iAM A 0T )Q k - Ql(G 00 + zAMA 0T )P + e mkn Ql G 00 Q n ] 



gl [q[ H 10 Q l - Q\ H 10 Q k + e klm (qI H 10 P - pt H 10 Q m )]} + oQ 



M 2 ") ' 
(4.11) 

where AM = Mq — M p and flavour indices have been suppressed. For AM = each row is 
separately invariant under the heavy quark spin symmetry, eq. (|4.2H . We focus on two of the 
operators in the following, namely ~ P^FP (cf. eq. (14.24p ) and ~ ig^P^GQ (cf. eq. (14.19P ). 

4.2. Euclidean formulation 

Before proceeding with the main line of the computation, we need to confront the specific time 
ordering appearing in eq. (|2.6j) . We do this by going through the imaginary-time formalism, 
because this produces an intermediate result which could in principle be compared with 
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lattice simulations. Denoting by r = it the Euclidean time coordinate, the propagator of P a 
is then determined by the Euclidean action (corresponding to a weight exp(— Se)) 

Sf = f dr f 2 Pj(«? T + M P )P a + (P a ~ Q ak ) . (4.12) 

JO Jx 

Due to charge conservation the susceptibility of eq. (|2.4j) can be re-expressed as 

r/9 



X 00 



J o dr J\[2Ptp a ](r,x) [2PjP 6 ](0,O)> T + (P a - Q ak ) + oQ^) , (4.13) 



and, taking into account the Wick rotation of the time coordinate, the Euclidean correlator 
related to eq. (|2.6f) can be defined as 

G E (r) = £ Hm -go <5°.7*(t> x) d°J k (0, 0)> T . (4.14) 

The subscript in G E could refer to "electric". In the subsequent steps we again omit hats, 
because eq. (|4.14p can be evaluated through normal Euclidean path integrals. 

The correlator of eq. (|4.14p is periodic across the Euclidean time interval, G E (r + k/3) = 
G e (t), k 6 Z, and can be Fourier analyzed in the Matsubara formalism. In particular, after 
a Fourier transform, 

G E (u n ) = f dre^G B (r) , (4.15) 
Jo 

where uj n = 2nnT, n 6 Wi, the spectral function is obtained from an imaginary part |23j : 

Pb(uj) = ImG E (w n -> -i[oj + i0 + ]) . (4.16) 



The momentum diffusion coefficient then follows from eq. (|2.7p . 

The Euclidean action in eq. (|4.12p implies that, for < |t — cr| < /3, the free propagator is 

<P a (T,x)P b W)>(°) 

e iu>n(T-cr) 



k^ 3 )(x-y)T£. 
\b ah ^X* - y)n B (M P )[e(a - r)e^ M r + 9{r - a)e^- T+ ^] , (4.17) 



where n B (M p ) = \j{eP M p — 1) is the Bose distribution. The propagator {Q a kQbi)T ^ as ^ e 
same structure, with an additional 5ki- The susceptibility of eq. (|4.13p then evaluates to 



X 00 = f}N t 5®(0) (e-^+3e-^)+0(^,^) , (4.18) 
where we also took the limit M p ,Mq » T, replacing n B (M) through exp(— (3M). 
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(a) (b) (c) 

Figure 2: Euclidean correlators corresponding to the scattering processes shown in fig. [2] The solid 
(single or double-lined) circle represents the Euclidean time interval; open squares correspond to force 
operators; and P and Q stand for the pseudoscalar and vector fields, respectively (in particle language 
the Q-field is denoted by P*). The pion lines do not meet in process (c). 



4.3. Leading order 

We now move on to consider the correlator G E , denned by eq. (|4.14p . with the force inserted 
from eq. (|4.1ip and the susceptibility from eq. (|4.18p . Some of the relevant Feynman graphs 
are shown in fig. [21 we start by considering the process (a). It was already argued in connection 
with fig. Q]that this should give no contribution to k but, as a preparation for the next-to- 
leading order computation, we recall briefly how the vanishing can be seen. 

Making use of eqs. (|3.8|) . (|4.8j) and the 2nd and 3rd rows of eq. (|4.11|) . the process in 
fig- corresponds to 

2 

(cPj k (r,x)cPj k (0,O)) T ~ j$pe([Ql A](r,x) [Pj ddfa Q d ](0, 0)) r , (4.19) 

where space-time indices have been omitted. Contracting the heavy-meson fields as in 
eq. (|4.17p , we note that there is one appearance of "forward" and one of "backward" propaga- 
tion; in the imaginary-time formalism, these amount to circling the Euclidean time interval. 
If we also insert the Goldstone boson propagator from eq. (13.2ft . the result becomes 



1 



(d°j k (T, x ) d°j k (o,o)) T 

- 2 .(iV f 2 
M 2 F 2 

where 



l l ,j(3) (o) n B (M P ) n B (M Q ) e ^ M P e M Q ddddA(r, 0) , (4.20) 



r „iu) n T+ip-0 r „ t rp \ r -, 



with E p = >\J m 2 + p 2 and J = Jd 3 p/(27r) 3 . A Fourier transform (cf. eq. (|4.15p ) can now be 
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taken: 



f 

Jo 



^ Te ^nr e (P-T)M Pe rM Q X T Ep ± e ^-r)E p 

e /3(M Q +E p ) _ e j3M P e PM Q _ e f)(M p +E v ) 



(4.22) 



iuj n + Mq — M p + E p iuj n + M Q - M p - E p 

The subsequent discontinuity, eq. (|4.16p . turns these into energy conservation constraints: 

e P(M Q +E p ) _ e j3M p e /3M Q _ e f3(M p +E p ) 



/j e (cl>) oc Im 
Jt 



p 1 uj + M q - M p + E p + iO+ ~ co + M Q - M p - E p + i0+ 



-7T 



| [<5(oj + M Q -M P + E p )e pEp ± 5(uj + M Q - M P - E p ) e m Q ( 



(4.23) 



where the ^-functions were made use of for rewriting the exponentials. (Using Hermitean 
conjugate operators in eq. (|4.19p yields the same but with Mq «-> M P .) 

It can now be seen that the limit of eq. (|2.7p can be non-trivial only if there are pion 
momenta with E p = +(Mq — M p ) (in fact this is a necessary but not a sufficient condition; 
in principle it could happen that the numerator of eq. (|4.20p . which has been left implicit here, 
vanishes at the same point, however this does not appear to be the case). This corresponds 
indeed to fig. [Ha), and cannot be realized if \Mq — M P \ < m n , as is the case with £>-mesons. 
Let us end by remarking that the same is expected to be the case also for many corrections 
of 0(g^/F 4 ), obtained from the topology of fig. E^a) by dressing either the pion propagator 
or the force operator by a closed "bubble" , which does not affect momentum flow. 

4.4. Next-to-leading order 

Moving on to next-to-leading order, we need to consider graphs like (b) and (c) in fig. [2l 
Due to the large number of terms proportional to in eq. ()4.1ip . as well as the possibly 
associated uncertainty as discussed below eq. (|4.2p . we simplify the task here by setting 
g w — > 0, omitting thereby contributions of O(g^) and O(g^) from k, which would otherwise 
be of the same order in the chiral expansion. Recalling that phenomenologically g n ~ 0.5, 
this cannot change the overall magnitude of the result. In this limit only the process (b) (as 
well as the same topology with P a — » Q ak ) is left over. Apart from "trivial" mass effects, spin 
symmetry implies that the terms containing P a and Qafc have identical structures at 0(g®). 
Therefore in the following only the part involving P a is displayed explicitly. 
As a first step we rewrite the relevant term of eq. (14. lip : 

d °[ jk = ^l P a( V >t-^l)n + 0(^) . (4.24) 
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This can be contrasted with eq. (2.17) of ref. [7J: the colour-electric field strength of HQET 
has been replaced by a kind of "chiral-electric" field strength in HM^PT. 
Going then over to the Euclidean correlator of eq. (|4.14p . we obtain 

G H (r) = | Jml ^({pl[vT,Vl] ab P b }(T,x){pM ,Vl] cd P d }(0,O)) T + 0(gl) . 

(4.25) 

Inserting the heavy meson propagator from eq. (|4.17j) and the susceptibility from eq. (|4.18p . 
this can be re-expressed as 

G E (r) = ^-([^,^] ab (r,0)[^,Pf] fea (0,0)) T + o(||,-^) . (4.26) 

With the contribution of Q ak added as has already been done here, the Boltzmann weights 
(cf. eq. (|4.18p ) have duly cancelled between the numerator and the denominator. 

Having obtained eq. (|4.26p . the problem has reduced to one within normal Chiral Pertur- 
bation Theory, eq. (f3~l~]h Noting from eqs. (pT8|) . (|4?7|) that 



[V^VZU = i{d T V k -d k V T ) ha + o(^) (4.27) 

= ^(d^d^-d^dkOba + oQ^, (4.28) 



and inserting the propagator from eq. (|3.2p . we obtain after some algebra that 

G E (r) = ^l[a T a fc A(r,0)5 T 5 fc A(r,0) - ^A(r, 0)V 2 A(r, 0)] +o(&,±r) . (4.29) 

Employing A(r, 0) from eq. (|4.2ip . the first term of eq. (|4,29p is seen not to contribute, 
because the p-integrand is odd in p — > —p. The second term of eq. (|4.29p does contribute; 
the Fourier transform (cf. eq. (I4.15|) ) can be carried out and, apart from a contact term oc S(r), 
yields a sum of four terms, with w n -dependence in structures of the type ~ l/(iu n + E p ±E q ), 
in analogy with eq. (14.221) . The spectral function (cf. eq. (|4.16p ) is obtained by replacing 
these with — tt5(u + E p + E q ), in analogy with eq. (I4.23j) . Thereby we end up with 

tt{N?-1) f p 2 E q 



J 

J p,q 



6F 4 L n 4E p 



x |[l + n B {E p ) + n B (E q )] [s(lo - E p - E q ) - 5(u + E p + E q ) 
+ [n B (E p ) - n B (E q )] + E p - E q ) - 5(u - E p + E q )} } 

+ 0(f|i). (4.30) 

The final step is to take the limit defined in eq. (|2.7p . Due to vanishing phase space, only 
the second structure of eq. (I4,30p gives a contribution linear in uj at u — » 0; this corresponds 
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to the process in fig. d{b) in the limit of vanishing energy transfer. The terms linear in u can 
be extracted by using the <5-functions to rewrite the arguments of the Bose distributions, and 
by then Taylor-expanding the latter. There are two terms, amounting to a symmetrization 
p <-» q; afterwards we can set E q — > E p . Given that n' B (E) = —f3n B (E)[l + n B (E)], this 
yields 



l±l nB(Ep) [l + n B (E p) } S (E p - E,) + ij) 



= ^ffV^^ + ^JM^ + O^,^), (4.31) 

where a partial integration was carried out. In the limit ttT » m w the integral can be 
performed explicitly and the result reads 



_ »T»m» 2(iV f 2 - l)vr 3 T 7 



l-^^+0(^A\+0( 9 ^ (4 32) 



63F 4 

A numerical evaluation is shown in fig. [3l This constitutes our final result. 



5. Conclusions and outlook 

The purpose of this paper has been to estimate the overall magnitude of the momentum 
diffusion coefficient, k, of a heavy-light pseudoscalar meson almost at rest with respect to a 
heat bath, which has a temperature in the range of some tens of MeV. The result, eq. (|4.3ip 
and fig. El is given in terms of low-energy constants of two-flavour Chiral Perturbation Theory. 

Although the result obtained is not immediately applicable to heavy ion collision experi- 
ments, in which the initial temperature may be in the range of hundreds of MeV and much 
of heavy quark energy loss could take place in the deconfined phase, the hope is that the 
analysis is nevertheless of theoretical interest. After all, physical QCD is believed to possess 
a crossover rather than a genuine phase transition between low and high temperatures, so 
that analyses in the former regime may yield qualitative information also for the latter. More 
concretely, our analysis combined with existing weak-coupling computations [6j suggest a 
picture in which the momentum diffusion coefficient exhibits a broad peak around the QCD 
crossover. So, heavy quark jets produced in an initial hard process may continue to approach 
kinetic equilibrium all the way until their final decoupling within the hadronic phase. 

Despite the qualitative nature of our study, it seems that a number of potentially interesting 
extensions can be envisaged. Most obviously, the result for k of 5-mesons could be completed 
with terms of 0(g%) and O(g^). Although no qualitative changes are expected, the corrections 
may be numerically important, perhaps in the range of ~ 50%. In addition the analysis 
appears formally interesting: as discussed in the text, it should probably be carried out both 
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Figure 3: An illustration of k/T 3 as a function of the temperature. The dashed line shows the 
result from eq. (|4.31|) : the dotted lines show two orders of the high-temperature expansion from 
eq. (|4.32l) . At high temperatures, T > 80 MeV, the chiral expansion breaks down because higher 
order corrections of relative magnitude 0(T 2 /F 2 ) become large; at low temperatures, T < 30 MeV, 
the high-temperature expansion breaks down because the expansion parameter 0{m 2 l /iT 2 T 2 ) exceeds 
unity (the unexpanded result of eq. (|4.31[) remains valid but is exponentially small). In the deconfincd 
phase a value k/T 3 > 2 might be phenomenologically acceptable (see e.g. ref. [15]; an oft-cited 
"diffusion coefficient" is D ~ 2T 2 /k < 1 /T) . The weak-coupling expansion suggests in general values 
k/T 3 < 2 [6], and at very high temperatures k/T 3 decreases like ~ l/ln 2 (T/To). 

within the non-relativistic theory as well as within a relativistic extension thereof, taking the 
non-relativistic limit only afterwards in the latter case, in order to have a crosscheck. 

Another interesting topic is the "extension" of the study to D-mesons. This immediately 
leads to the dramatic effect that resonant contributions, such as the process shown in fig.QJa), 
are allowed, whereby the .D-mesons can change their identity (a further complication is that 
for D* decays electromagnetic processes are important). The handling of the related rich 
physics may suggest an excursion away from systematic computations, into the realm of 
hadronic models; an overview of recent works in this direction can be found e.g. in ref. |40| . 
Nevertheless, it can be noted that resonant phenomena could also play a role in theoretical 
considerations of the U-sector, if we took the chiral limit — > while keeping the mass 
difference M Bn — M B fixed and non-zero. 

Yet another line is that whereas only the "intercept", n = lim^^o 2Tp E (a;)/u;, was addressed 
here, the full spectral function p E (^) (eq. (14. 16ft ). and even the Euclidean correlator G E (r) 
(eq. (|4.14p ). are also interesting objects. For instance, they may help in the interpretation 
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of the corresponding lattice measurements a la ref. [19] (if these were unquenched); it is in 
this spirit that both functions have been computed also in the deconfined phase [18], and in 
fact G e (t) even in the confined phase in the presence of a large lattice spacing, through the 
use of the lattice strong-coupling expansion [18] . The terms of 0(g^/F 2 ) (cf. fig. [2T a )) do 
contribute significantly to G E (r), even if they do not contribute to k as has been discussed in 
the text; it might turn out to be useful to gain understanding on how hidden the information 
about k is in the directly measurable G E (r). 

A further topic is not to consider the spectral function p E (w) yielding k, but rather the full 
spectral function /o(w,k) = L x e 4 "^*'*^ [<^ M (*> x )i ^t(O) 0)]) T , of which p E (w) is a specific 
limit [7j. Restricting first to vanishing momentum, k = 0, and to small frequencies, the 
analysis can still be carried out within the non-relativistic framework. This spectral function 
is expected to show a narrow and prominent "transport peak" at frequencies \cj\ < rj D ~ 
k/(2TM). A direct analysis of such infrared features tends to be difficult, necessitating 
complicated resummations (see e.g. ref. [22])) but at least the theory in question is not a 
gauge theory, whereby these may be more tractable than in the deconfined phase. 

The most ambitious goal would be to consider p(ui, k) also for k # 0, perhaps even for 
|k| > M, corresponding to heavy mesons moving at a relativistic speed with respect to a 
pionic plasma. Although this could take us away from the range of validity of systematic 
HM^PT, simple relativistic extensions can be written down (see e.g. refs. [37, 39J). Within 
such a framework one might also try to understand phenomena such as radiative energy loss, 
perhaps making contact with classic pion gas computations that have recently been revived 
through the AdS/CFT setup [41]. 
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Note added 

After the submission of this paper three works appeared |42]-|44] in which the kinetic equi- 
libration of D-mesons is studied in a spirit similar to that suggested in sec. [5j 
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